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■ ABSTRACT 

I : Context. The study of the younger, and brighter, pulsars is important to understand the optical emission properties of isolated 
O ' neutron stars through observations which, even in the lOm-class telescope era, are much more challenging for older and fainter 
£j objects. PSR B0540— 69, the second brightest (V ~ 22) optical pulsar, is obviously a very interesting target for these investi- 
ffi ■ gations. 

i Aims. The aim of this work is threefold: constraining the pulsar proper motion and its velocity on the plane of the sky, obtaining 

a more precise characterisation of the pulsar optical spectral energy distribution (SED) through a consistent set of multi-band, 
y—i • high-resolution, imaging photometry observations measuring the pulsar optical phase-averaged linear polarisation, for which 
^ ' only a preliminary and uncertain measurement was obtained so far from ground-based observations . 

Methods. We performed high-resolution observations of PSR B0540— 69 with the Wide Field and Planetary Camera 2 (WFPC2) 
aboard the Hubble Space Telescope (HST), in both direct imaging and polarimetry modes. 

Results. From multi-epoch astrometry we set a 3a upper limit of 1 mas yr _1 on the pulsar proper motion, implying a transverse 
velocity < 250 km s -1 at the 50 kpc LMC distance. Moreover, we determined the pulsar absolute position with an unprece- 
C*") ' dented accuracy of 70 mas. From multi-band photometry we characterised the pulsar power-law spectrum and we derived the 
most accurate measurement of the spectral index (qo = 0.70 ± 0.07) which indicates a spectral turnover between the optical 
and X-ray bands. Finally, from polarimetry we obtained a new measurement of the pulsar phase-averaged polarisation degree 
I (PD = 16% ± 4%), consistent with magnetosphere models depending on the actual intrinsic polarisation degree and depolari- 
^» , sation factor, and we found that the polarisation vector (22° ± 12° position angle) is possibly aligned with the semi-major axis 
! *~j ■ of the pulsar- wind nebula and with the apparent proper motion direction of its bright emission knot. 

r> | Conclusions. By using the WFPC2 on the HST, we performed a comprehensive optical study (astrometry, photometry, and 
^ polarimetry) of PSR B0540— 69. Deeper studies with the HST can only be possible with the refurbished Advanced Camera for 
■ ■ • Surveys {ACS) and with the new Wide Field Camera 3 (WFCZ). 
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1. Introduction 

The pulsar PSR B0540-69 in the Large Magellanic Cloud 
(LMC) is often referred to as the Crab pulsar "twin" be- 
cause of it is very similar in age (~ 1700 years), spin pe- 
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riod (P = 50 ms), and rotational energy loss (E ~ 10 38 
erg s _1 ). PSR B0540— 69 is the second pulsar discovered 
in X-rays by the Einstein Observatory (Seward et al. 1984) 
and the first extragalactic pulsar detected at any wave- 
length. Like the Crab, PSR B0540-69 is also embedded in 
a bright pulsar wind nebula (PWN) visible at wavelengths 
from the optical to the soft/hard X-rays (e.g. De Luca et 
al. 2007; Petre et al. 2007; Slowikowska et al. 2007). After 
its discovery, PSR B0540— 69 has been observed by nearly 
all X-ray satellites, EXOSAT (Ogelman & Hasinger 1990), 
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Ginga (Nagase et al. 1990), ROSAT (Finley et al. 1993), 
BeppoSax (Mineo et al. 1999), Chandra (Kaaret et al. 
2001), Rossi-XTE (de Plaa et al. 2003), ASCA (Hirayama 
et al. 2002), Integral (Slowikowska et al. 2007), and Swift 
(Campana et al. 2008). Its X-ray light curve is charac- 
terised by a single, broad peak, very much at variance 
with, e.g. that of the Crab pulsar. In radio, the distance 
to the LMC made PSR B0540-69 undetectable for a long 
time until pulsations were finally detected from Parkes 
(Manchester et al. 1993). Giant radio pulses were discov- 
ered by Johnston & Romani (2003), aligned in phase with 
the peak of the X-ray pulse (Johnston et al. 2004), making 
PSR B0540— 69 the second youngest radio pulsar to fea- 
ture this phenomenon. PSR B0540— 69 is also one of the 
handful of rotation-powered pulsars with a measured brak- 
ing index, obtained from X-ray observations (e.g., Zhang 
et al. 2001; Cusumano et al. 2003; Livingstone et al. 2005). 

In the optical, PSR B0540-69 is the second bright- 
est isolated neutron star (V ~ 22) identified so far (see 
Mignani 2009a, b for recent reviews). Optical pulsations 
were detected by Middlcditch & Pcnnypacker (1985) soon 
after the X-ray discovery, making PSR B0540— 69 the 
third optical pulsar after the Crab (Cocke et al. 1969) 
and Vela pulsars (Wallace et al. 1977). However, it was 
only through high-resolution imaging observations with 
the ESO New Technology Telescope (NTT) that its op- 
tical counterpart was actually identified (Caraveo et al. 
1992; Shearer et al. 1994). The PSR B0540-69 optical 
light curve (Middleditch et al. 1987; Gouiffes et al. 1992; 
Boyd et al. 1995; Gradari et al. 2009) is characterized by 
a single broad peak, very similar to the X-ray and ra- 
dio light curves, with a significant dip on top. The optical 
spectral energy distribution (SED) of PSR B0540-69was 
first measured by Middleditch et al. (1987) from high- 
speed multi-band photometry which suggested the pres- 
ence of a possible excess in the U band with respect to an 
otherwise monotonic power-law continuum (F v oc v~ a °). 
A re-analysis of the same multi-band photometry mea- 
surements, however, did not yield evidence of the claimed 
U-band excess (Nasuti et al. 1997). The power-law spec- 
trum demonstrates that the optical emission from PSR 
B0540— 69 is of magnetospheric origin, with an optical lu- 
minosity consistent with the expectations of the Pacini 
& Salvati (1987) model. Low-resolution spectra of PSR 
B0540-69 were obtained by Hill et al. (1997), with the 
Hubble Space Telescope (HST), and by Serafimovich et 
al. (2004), with the Very Large Telescope (VLT), although 
the latter was affected by the contamination from the su- 
pernova remnant. By using multi-band imaging photom- 
etry of PSR B0540— 69 from archival HST observations, 
Serafimovich et al. (2004) confirmed that the optical spec- 
trum is dominated by a power-law continuum, although 
with a spectral index different from the values previously 
published. They also reported a tentative proper motion 
measurement for PSR B0540— 69, which, however, was not 
confirmed by De Luca et al. (2007) using a large time 
baseline of HST observations. Phase-resolved polarime- 
try observations of PSR B0540— 69 were performed by 



Middleditch et al. (1987) but only yielded an upper limit 
on the phase-averaged polarisation degree, while image 
polarimetry observations of Chanan & Helfand (1990) fo- 
cused on the PWN only. More recently, image polarime- 
try observations of PSR B0540— 69 were performed by 
Wagner & Seifert (2000) using the VLT. They reported 
a phase-averaged polarisation degree of ~ 5% (with no 
associated error) which was probably affected by the con- 
tribution of the unresolved PWN. Thus, only preliminary, 
or uncertain, phase-averaged optical polarisation measure- 
ments exist for PSR B0540-69. 

In this paper, we report the results of new astrome- 
try, photometry, and polarimetry measurements of PSR 
B0540— 69 performed with the HST as a part of a dedi- 
cated programme aimed at the comprehensive study of the 
pulsar and of its PWN (Mignani et al., in preparation). 
From the pulsar astrometry we will constrain its proper 
motion and its velocity on the plane of the sky, and we 
will provide an update reference position for future ob- 
servations. From the pulsar photometry we will measure 
anew the optical SED, we will study the relation with 
the X-ray spectrum, and we will verify the presence of 
a spectral turnover between the two energy bands. From 
the pulsar polarimetry we will measure its phase-averaged 
polarisation and we will test the predictions of different 
neutron star magnetosphere models. The paper is divided 
as follows: observations and data analysis are described in 
Sect. 2, while the results are presented and discussed in 
Sect. 3 and 4, respectively. 

2. Observations and data analysis 

We observed PSR B0540-69 with the Wide Field and 
Planetary Camera 2 (WFPC2) aboard HST (Programme 
#10900, PI Mignani). The WFPC2 is a four-chip CCD 
detector, sensitive to radiation in the 1150 — 11000 A spec- 
tral range. In order to exploit the maximum spatial resolu- 
tion, the pulsar was centred in the Planetary Camera (PC) 
chip which has a pixel scale of 0'.'045 and a field-of-view 
of 35" x 35". The WFPC2 observations were performed 
on June 21st, July 26th, September 25th, and November 
5th, 2007 for a total of six spacecraft orbits, correspond- 
ing to four different visits. We note that since our target is 
in the continuous viewing zone (CVZ), a declination zone 
which can be visible for an entire spacecraft orbit for a few 
days of the year, our programme could have been carried 
out in one visit only. However, we decided to split it into 
different visits since the narrow CVZ time window would 
have made more difficult to cope with the spacecraft ori- 
entation constraints required for the WFPC2 polarimetry 
observations (see below). 

As part of our programme, we performed both multi- 
band photometry and polarimetry observations (see Table 
[TJ top section, for a summary). To maximise the spec- 
tral coverage, we observed through the F336W, F450W, 
F555W, F675W, and F814W filters. Observations were 
carried out in the first visit for two consecutive orbits of 
the spacecraft. The last orbit of the first visit was devoted 
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Table 1. Summary of the HST/WFPC2 observations 
PSR B0540-69, with the observing dates, the filter 
names, their pivot wavelength A and bandwidth A A in 
A, and the total integration time T in seconds. 



Date 


Filter 


A (AA) 


T 


yyyy-mm-dd 




(A) 


(fl) 


2007-06-21 


F336W 


3359 (204) 


780 




F450W 


4557 (404) 


780 




F555W 


5442 (522) 


300 




F675W 


6717 (368) 


420 




F814W 


7995 (646) 


600 




F606W/POLQ 


6038 (651) 


1800 


2007-07-26 


F606W/POLQ 




1800 


2007-09-25 


F606W/POLQ 




1800 


2007-11-05 


F606W/POLQ 




1800 


2005-11-15 


F555W 




480 




F547M 


5483 (205) 


1040 


1999-10-17 


F336W(*) 




600 




F791W(*) 


7872 (519) 


400 




F547M(*) 




800 




F502N 


5012 (27) 


11000 




F673N 


6732 (47) 


8200 


1995-10-19 


F555W(*) 




600 




F656N 


6564 (22) 


3400 




F658N(*) 


6591 (29) 


4000 



(*) observations used by Serafimovich et al. (2004) 



to the first of the four planned polarimetry observations, 
one at each of the selected polarimetry angle. Indeed, the 
only possibility to perform polarimetry observations with 
the WFPC2 and to keep the target positioned in the same 
chip, is to rotate the spacecraft to orient the instrument 
polariser. This requires that, for each polarimetry angle, 
the observations must be scheduled in different visits. This 
means that the choice of the spacecraft roll angle is lim- 
ited by the guide star/solar panel constraints. Polarimetry 
observations were performed at nominal position angles of 
0°, 45°, 90°, and 135° using the F606W filter in combina- 
tion with the Polariser Quad (POLQ) filter. For each filter, 
all observations were split into three shorter exposures to 
filter out cosmic ray hits. 

Additional WFPC2 observations of PSR B0540-69, 
taken through broad-band filters, are available in the 
HST archive (see Table [TJ lower sections) . These include 
two sets of F555W observations taken on October 19th 
1995 (#6120, PP Caraveo) and November 15th 2005 
(#10601, PP Lundqvist), respectively, and one set of 
F336W and F791W observations, taken on October 17th 
1999 (#7340, PI: Morse). Observations taken through 
the medium-band filter F547M are available from pro- 
grammes #7340 and #10601. Finally, observations taken 
through the narrow-band filters F656N, F658N, F502N, 
and F673N are available from programmes #6120 and 
#7340. Part of the 1995-1999 WFPC2 observations were 
used by Serafimovich et al. (2004) to study the pulsar op- 
tical SED. To complete the available spectral coverage, we 



decided to add to our data set the archival observations 
taken through the F791W filter and, for comparison, those 
taken through the F336W and F555W ones. At variance 
with Serafimovich et al. (2004), we decided not to use the 
observations taken through narrow-band filters, where the 
pulsar is fainter and where it is more difficult to isolate its 
emission from that of the nebula. 

All data were processed on-the-fly at the Space 
Telescope Science Institute (STScI) through the WFPC2 
CALWP2 reduction pipeline for bias, dark, flat-field correc- 
tion, flux, and polarimetry calibration using the closest in 
time reference calibration frames. For each filter, we com- 
bined single exposures with the STSDAS task combine to 
produce co-added and cosmic-ray free images. 

3. Results 

3.1. Astrometry 

Firstly, we used the available high-resolution WFPC2 ob- 
servations of PSR B0540— 69 to obtain a compelling con- 
strain on the pulsar proper motion through optical as- 
trometry. To do this, we used the long time baseline pro- 
vided by the observations of De Luca et al. (2007). In par- 
ticular, we used all data collected with the F555W and 
F547M filters (see Table [IJ and we considered only the 
PC images to maximise the spatial resolution. The se- 
lected data set includes five observations spanning ~ 11.7 
years. In order to measure the PSR B0540— 69 proper mo- 
tion through relative astrometry, we applied the algorithm 
that we have successfully applied in a number of previous 
investigations including, e.g. our first study of the pul- 
sar proper motion (De Luca et al. 2007). We started by 
selecting a grid of 60 good reference sources which were 
detected in the PC ficld-of-view and not extended, not 
saturated, with a high signal-to-noise ratio and not too 
close to the CCD edges. For each source, we computed 
their pixel coordinates on each image by fitting a 2-D 
gaussian to their brightness profile, with an uncertainty 
of 0.02-0.06 pixels per coordinate. We evaluated the po- 
sition of the pulsar optical counterpart in the same way, 
with an uncertainty of 0.03-0.04 pixels per coordinate. We 
corrected the source coordinates for the WFPC2 geomet- 
ric distortion (Anderson & King 2003), as well as for the 
u 3A th row defect" (Anderson & King 1999). We assumed 
the 1995 image as a reference and we aligned the associ- 
ated reference grid along right ascension and declination 
using the known telescope roll angles. By fitting the refer- 
ence star positions, we computed the coordinate transfor- 
mation which yields the best superposition of each frame 
grid on the 1995 one. In the fitting procedure we discarded 
ten reference object which yielded large residuals using an 
iterative sigma clipping algorithm. The resulting r.m.s. un- 
certainty on the superposition of the frame grids turned 
out to be 0.06-0.08 pixel per coordinate (using the remain- 
ing 50 reference stars). This coordinate transformation 
allowed us to convert the pulsar positions to a common 
reference frame and to evaluate its displacement in right 
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ascension and declination with respect to its 1995 posi- 
tion. A simple linear fit to the measured displacements 
as a function of time yields no evidence for any signifi- 
cant proper motion of the pulsar, with a 3er upper limit of 
0.015 pixel yr _1 per coordinate. Using the well calibrated 
WFPC2 plate scale, we set a 3a upper limit to the overall 
pulsar proper motion of 1 mas yr _1 . Such a result con- 
firms our previous findings (De Luca et al. 2007), but it is 
slightly more compelling because of the longer time base- 
line (11.7 years with respect to ~ 10 years) covered by the 
used data set. 

Secondly, we used our new WFPC2 images to com- 
pute updated coordinates for PSR B0540— 69 through 
optical astrometry. To this aim, we re-calibrated the im- 
age astrometry using as a reference the position of stars 
selected from the Two Micron All Sky Survey (2MASS) 
catalogue (Skrutskie et al. 2006) which has a mean po- 
sitional error Ss 0'.'2 and it is linked to the International 
Celestial Reference Frame (ICRF) with a ~ 0'.'015 accu- 
racy (Skrutskie et al. 2006). As discussed in Mignani et al. 
(2005), 2MASS is preferred to the Guide Star Catalogue 
2.3 (GSC-2.3; Lasket er at. 2008) and to the U.S. National 
Observatory B1.0 catalogue (USNO-B1.0; Monet et al. 
2003) for fields in the LMC. On the other hand, both the 
USNO CCD Astrograph Catalogue 2 (UCAC-2; Zacharias 
et al. 2004) and its update UCAC-3 (Zacharias et al. 
2009) only provide a sparse astrometric grid in the mosaic 
WFPC2 field-of-view. As a reference for our astrometry re- 
calibration, we chose the June 2007 image taken through 
the F555W filter, where PSR B0540-69 is detected with a 
good signal-to-noise, and for which the correction for the 
WFPC2 geometric distortion is well modelled. We used 
the mosaic of the four WFPC2 chips (160" x 160" field- 
of-view; O'.'l pixel scale) since it provides a large enough 
ficld-of-view to include a sufficient number of 2MASS stars 
to be used for the astrometry re-calibration. We pro- 
duced the mosaic image with the STSDAS task wmosaic 
which also applies the correction for the geometric dis- 
tortions of the four chips. Approximately 60 2MASS ob- 
jects are identified in the mosaic WFPC2 image. As done 
in Sect 3.1, we filtered out extended objects, stars that 
are either saturated or too faint to be used as reliable 
astrometric calibrators, or too close to the CCD edges. 
We finally performed our astrometric calibration using 40 
suitable 2MASS reference stars, evenly distributed in the 
WFPC2 field-of-view. Then, for both the 2MASS stars 
and the pulsar, we determined their pixel coordinates in 
the WFPC2 reference frame from the centroid of their 
intensity profile, computed by fitting a 2-D Gaussian. 
This yielded errors (per coordinate) of Si 0.01 and of 
< 0.05 WFPC2 pixels on the centroids of the 2MASS 
stars and of the pulsar, respectively. We then computed 
the fit to the pixel-to-sky coordinate transformation for 
the 2MASS stars using the code ASTROrfl based on higher 
order polynomials. By applying the computed fit to the 
pixel coordinates of PSR B0540— 69 we finally obtained: 

1 http:/ /star-www.rl.ac.uk/Software/software.htm 



aj2ooo = 05 h 40 m 11.204 s , Sj 20 oa = -69°19'54"34, with an 
overall accuracy Sr — 0'.'12 (la). We estimated the accu- 
racy on the computed coordinates by adding in quadrature 
the rms of the astrometric fit, ar ~ O'.'l, the uncertainty 
in the registration of the WFPC2 image on the 2MASS 
reference frame, a tr = 0"05dl, the ~ 0"015 accuracy of 
the link of the 2MASS coordinates to the ICRF, and the 
uncertainty on the pulsar centroid (Ss 0'.'005). Since the 
uncertainty on the reference star centroids is much smaller 
than that of the pulsar, we neglected it in our global error 
budget. 

As a check of our absolute astrometry, we used each 
of the F555W data sets (see Table [I}, which were ob- 
tained from observations performed in different visits, i.e. 
with different telescope pointings and spacecraft roll an- 
gles, to obtain independent measurements of the pulsar 
position. Since the pivot wavelength of the F547M filter is 
essentially the same as that of the F555W one (see Table 
[I]), the wavelength dependence of the geometric distor- 
tion should not induce any bias when using the correc- 
tion optimized for the F555W filter (see also De Luca et 
al. 2007). Thus, we could also use the two F547M data 
sets as a reference for our absolute astrometry. In par- 
ticular, for the November 2005 observations, we decided 
to use the F547M data set instead of the F555W one, 
due to its longer integration time. We note that the other 
suitable medium/broad-band filter data sets can not be 
formally used for an independent check of our astrome- 
try since they were obtained from observations performed 
during the same visits as the F555W and F547M ones and, 
thus, with the same telescope pointing and roll angles. In 
principle, the data sets obtained from the observations 
performed through the F606W/POLQ filter, which were 
taken in different visits, could be used, though. However, 
possible effects of the polarisation optics on the WFPC2 
astrometry, as well as the mapping of the geometric dis- 
tortion, are still to be studied, thus making the use of 
imaging polarimetry data less suitable for astrometry (see 
also Kaplan et al. 2008). By re-calibrating the astrom- 
etry of each data set through the same procedure de- 
scribed above we obtained: aj2ooo = 05 ft 40 m l 1.208 s , and 
<Sj20oo = -69°19'54"12 (Sr = 0"14; la) for the October 
1995 F555W data set. Similarly, we obtained: aj2ooo = 
05 /l 40 m l 1.193 s , <5j 20 oo = -69°19'54"12 (Sr = 0"15; la) 
and ajsooo = 05 /l 40 m 11.206 s , <5 J2 ooo = -69°19'54"10 
(Sr = 0"13; la) for the November 2005 and for the 
October 1999 547M data sets, respectively. We note that 
all the four sets of coordinates are consistent within s» la, 
which confirms that our first measurement was free of 
systematics. The marginal, although not significant, dif- 
ferences in the pulsar coordinates are most likely due to 
the different number of 2MASS stars used to compute the 



2 As in Lattanzi et al (1997), we defined atr = \/3x a s /yfN s , 
where v3 accounts for the free parameters in the astrometric 
fit (x-scale, y-scale, and rotation angle) N a is the number of 
2MASS stars and <j s = 0'.'2 is the conservative mean positional 
error of their coordinates (Skrutskie et al. 2006) 
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astrometric solution and to their relative distribution in 
the WFPC2 field-of-view, which depends on the different 
telescope pointing directions and roll angles. Marginal dif- 
ferences can also be due to time-dependent shifts of the 
WFPC2 chips in the HST focal plane, which can affect 
astrometry of the WFPC2 mosaics. We note that the up- 
per limit that we obtained on the PSR B0540— 69 proper 
motion (see Sect. 3.1) is much smaller than the error on 
our absolute astrometry. Thus, we can neglect any dis- 
placement of the pulsar in the 11.7 year time span cov- 
ered by the available WFPC2 observations and we can 
compute an even more precise position from the average 
of the four sets of coordinates independently computed 
from the F555W/F547M data sets. This yields aj 20 oo = 
05 ft 40 m 11.202 s , 5/2000 = -69°19'54"17 {Sr = 0'.'07; la) as 
the most accurate value of the PSR B0540— 69 position. 



3.2. Multi-band photometry 

We used our WFPC2 observations of PSR B0540-69 to 
perform multi-band photometry using, for the first time, 
a set of observations consistently taken at the same epoch 
(see Table [T]). This is important to minimise systematic 
effects on the pulsar photometry due, e.g. to possible long- 
term variations in the instrument zero points, to a de- 
grade in the instrument efficiency, or to the use of differ- 
ent calibration data sets. In the case of PSR B0540— 69, 
high-resolution imaging photometry observations are, at 
present, better suited than spectroscopy since they are 
less affected by the contamination from the bright super- 
nova remnant (see also discussion in Serafimovich et al. 
2004). 

We measured the flux of PSR B0540-69 in our 
WFPC2 images through customised aperture photom- 
etry using the IRAF package digiphot. To maximise 
the signal-to-noise ratio, we measured the pulsar counts 
through small apertures (3 pixel radius) and we subtracted 
the PWN background measured in an annulus of 10 pixel 
radius at a 2 pixel distance from the photometry aper- 
ture not to include the wings of the point spread func- 
tion (PSF) . We then applied the aperture correction using, 
per each filter, the specific coefficients given in Holtzman 
et al. (1995). Observations at shorter wavelengths with 
the WFPC2 are affected by the presence of contaminants 
which progressively build up on the CCD surface (Bagget 
et al. 2002). The CCD contamination results in tempo- 
ral variations in throughput which cause a decrease in the 
measured source flux. This is an issue for the measurement 
of the pulsar flux in the F336W filter, where the contam- 
ination rate can be of the order of ~ 0.06% per day and 
tends to be higher at the centre of each PC chip, where 
the pulsar is positioned. This effect is not accounted for 
by the CALWP2 pipeline and has to be corrected manually. 
We then applied to our measurement the contamination 
correction, computed during monthly de-contamination 
procedures of the CCD, the last one carried out about 



three weeks before our observation^. To all our flux mea- 
surements we then applied correction to compensate for 
the time and position-dependent charge transfer efficiency 
(CTE) losses of the WFPC2 detector, using the formulae 
given in Dolphin (2009). 

Following the WFPC2 Data Handbook (Bagget et al. 
2002), we computed the count-rate to flux conversion us- 
ing the image keywords PHDTFLAM and PH0TZPT, respec- 
tively, derived by the WFPC2 photometry calibration 
pipeline. We then derived magnitudes rriF336W = 21.77 ± 
0.17, to F 4 5 ow = 22.09 ± 0.05, m F555 w = 22.07 ± 0.05, 
™-F675W = 22.18 ± 0.04, and m F81iw = 22.35 ± 0.03. 
To complement our broad-band imaging photometry mea- 
surements, we also computed the pulsar flux from the 
October 1999 F791W observation (see Table P). Using the 
same approach described above, we derived niF79iw 
'22.2'.) ± 0.04, which is in quite good agreement with that 
obtained through the slightly redder F814W filter. The 
computed magnitudes, together with their errors, are sum- 
marised in Table [2] The quoted formal errors are purely 
statistical and do not account for systematic errors related 
to the uncertainty of the WFPC2 photometry calibration, 
which are of the same order of magnitude for all filters. 

Table 2. Observed WFPC2 magnitudes of PSR 
B0540— 69 and associated errors (in parenthesis). 
Magnitudes measured in the present work are listed in the 
third column while those measured by Serafimovich et al. 
(2004) are listed in the fourth column for comparison. For 
the latter, values accounting for the CTE correction are 
given in brackets. 



Filter 


Date 

yyyy-mm-dd 


Ma 


gnitudes 




F336W 


2007-06-21 


21.77 (0.17) 








1999-10-17 


21.96 (0.19) 


22.25 [21.98] 


(0.13) 


F450W 


2007-06-21 


22.09 (0.05) 






F555W 


2007-06-21 


22.07 (0.05) 








2005-11-15 


22.11 (0.04) 








1995-10-19 


21.99 (0.03) 


22.17 [22.14] 




F675W 


2007-06-21 


22.18 (0.04) 






F719W 


1999-10-17 


22.29 (0.04) 


22.44 [22.37] 


(0-12) 


F814W 


2007-06-21 


22.35 (0.03) 







To have an independent cross-check of our measure- 
ments, we also measured the pulsar flux from the archival 
October 1999 F336W observation and from the November 
2005 and October 1995 F555W observations. Not to intro- 
duce systematic effects in our photometry, for both filters 
we used the same apertures and background areas used 
for the June 2007 observations. For the F336W filter we 
obtained mp 336 w = 21.96 ± 0.19, i.e. very well consis- 
tent with our measurement from the June 2007 obser- 



http: / / www.stsci.edu/hst /wfpc2 / analys 
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vation. For the the F555W filter, however, we obtained 
m F 555W = 22.11 ± 0.04 and m F555W = 21.99 ± 0.03 
for the November 2005 and for the October 1995 ob- 
servations, respectively. Thus, while the June 2007 and 
November 2005 F555W measurements are in almost per- 
fect agreement with each other, the October 1995 one 
shows a somewhat larger difference, between ~ —0.08 and 
~ —0.12 magnitudes. This difference, not formally signif- 
icant though, is larger than the nominal uncertainty of 
0.02-0.04 magnitues on the WFPC2 zero points (Heyer 
et al. 2004) and, thus, it implies an actual difference in 
the measured CTE-corrected count-rate, possibly related 
to the on-the-fly re-calibration of the WFPC2 October 
1995 F555W data set. We also compared our measure- 
ments with the pulsar fluxes measured by Serafimovich 
et al. (2004) on the very same data sets (fourth column 
in Table 0), i.e. for the October 1995 F555W and for the 
October 1999 F336W and F791W observations. We note 
that our measurements yield fluxes brighter by ~ 0.1 — 0.2 
magnitudes with respect to those of Serafimovich et al. 
(2004), although still compatible within the errors. We 
attribute this difference to the fact that they apparently 
did not apply the CTE correction to their measurements, 
which can be of this order of magnitude. Indeed, by apply- 
ing such correction we indeed found that their measure- 
ments turn out to be more consistent with ours (see Table 
[21 values in brackets). In particular, their CTE-corrected 
October 1995 F555W flux measurement agrees with our 
June 2007/November 2005 F555W ones, thus confirming 
that the slightly deviant measurement that we obtained 
from the former data set might actually reflect a difference 
in the on-the-fly re-calibration. 

3.3. Polarimetry 

We measured the pulsar polarisation using the web inter- 
fac^_] to the WFPC2 polarisation calibration tool (Biretta 
& Mc Master 1997). This program computes the trans- 
mission of a polariser element, for a given filter and CCD 
gain, using a synthetic spectrum simulated by the STSDAS 
package synphot. As a reference for the pulsar, we as- 
sumed its power-law spectrum determined in Sect. 3.2. 
Then, the program computes the Stokes /, Q, and U pa- 
rameters from the counts measured in the three images 
oriented along the three different polarisation angles and 
applies the coordinate transformation to convert the mea- 
sured polarisation from the detector to the sky reference 
frame using the telescope roll angle information. 

We selected the triplet of images oriented along the 
three position angles of 0°, 45°, and 90°. As done in 
Sect. 3.2, for each image we measured the pulsar counts 
with the IRAF digiphot package using an aperture of 
3 pixel radius and we subtracted the background sam- 
pled in an annulus of 10 pixel width. Since the WFPC2 
calibration tool returns the Stokes parameters as a lin- 
ear combination of the counts measured in the three im- 

4 www.stsci.edu/hst/wfpc2/software/wfpc2_pol_calib.html 



ages (Biretta & Mc Master 1997), this procedures accounts 
for the subtraction of the PWN contribution (Chanan & 
Helfand 1990) to the measured pulsar polarisation. We 
then applied the aperture correction from Holtzman et al. 
(1995) and the CTE correction from Dolphin (2009) to 
the background-subtracted source counts. In this way, we 
measured the Stokes parameters of the pulsar {Q / 1) psr = 
0.09%±0.03%and (U/I) PSR = 0.12%±0.02%. These cor- 
responds to a phase-averaged polarisation degree PD = 
15.2% ± 2.5% with the polarisation vector oriented along 
a positional angle PA = 21.1° ± 10° (east from north), 
where PA and PD are defined as: 

PA=[{Q/If + {U/I) 2 ] 1 ' 2 (1) 

PD = i arctan(LVQ) (2) 

The associated uncertainties on PD and PA are de- 
rived from the statistical errors on the pulsar counts mea- 
sured in each of the three images. To investigate a possi- 
ble dependence of the polarisation degree on the telescope 
roll angle chosen to align the polariser (see Sect. 2), we 
repeated our measurement using a different triplet of im- 
ages taken with a different sequence of roll angles (45°, 
90°, and 135°) but we found no significant difference. To 
the statistical uncertainty of our measurement we have 
to add the uncertainty on the absolute polarimetry cali- 
bration of the WFPC2, which includes the effects of the 
instrument intrinsic polarisation. For the PC chip, this un- 
certainty is estimated to be of the order of 1.5% (Biretta 
& Mc Master 1997). We thus assumed an overall error of 
3% on the measured polarisation degree of the pulsar. 

Of course, this value must be corrected for the con- 
tribution of the foreground polarisation towards PSR 
B0540— 69, which is produced by the dust in the inter- 
stellar medium (ISM) and by the integrated components 
of both the Galaxy and the LMC and it is not negli- 
gible, in principle. To estimate the contribution of the 
foreground polarisation, we measured the polarisation de- 
gree of a number of test stars uniformly distributed across 
the PC field -of-view. We followed the same recipe used 
for the pulsar but through an automatic procedure for 
source detection and count measurement in each of the 
three images, aperture and CTE corrections, source list 
matching, and Stokes parameter computation through the 
WFPC2 polarisation calibration tool. For the test stars we 
assumed a flat spectrum to compute the transmission of 
the POLQ polariser, which is accurate within ~ 1% as 
first approximation (Biretta & Mc Master 1997). From 
the starting sample of test stars we then selected those 
which are brighter than mF555W = 22 and which can thus 
be used as a better reference for our estimate. From these 
stars we then selected those (~ 50) for which the polar- 
isation degree PD is measured with a significance of at 
least 3ct. After applying a ler clipping filter, we ended up 
with a working sample of ~ 20 test stars from which we de- 
rived the average Stoke parameters (Q/I) s tars = —0.029% 
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(0.023% rms) and {U/I) stars = 0.015% (0.029% rms). 
From these values we computed an average polarisation 
degree PD star s = 3.3% ± 2.5% and an average a position 
angle PA star s ~ 75°. After accounting for the ~ 1.5% 
uncertainty on the absolute polarimetry calibration of the 
WFPC2, we then ended up with PD stars = 3.5% ± 3%, 
which we assumed as an estimate of the foreground polar- 
isation. We note that our value is compatible with the 
w 2 - 3% estimated by Chanan & Helfand (1990) us- 
ing a larger sample of test stars. After subtracting the 
foreground polarisation in the Stokes parameter space, 
the corrected Stokes parameters of the pulsar are then 
(Q/I)psr — 0.12% ± 0.03% and (JJ/I)psr — 0.11% ± 
0.03%. From these, we finally computed the intrinsic po- 
larisation degree of the pulsar PD = 16% ± 4% and the 
polarisation position angle PA = 22° ± 12°, where the 
attached errors also account for the rms on the Stokes 
parameters of the stars used to estimate the foreground 
polarisation. 



Table 3. Compilation of the PSR B0540-69 coordinates 
available from the literature and associated errors Aa and 
AS. 



(htns) 
a J2000 


Aa (s » 


°J2000 


A5^ 


Source 


05 40 11.03 


0.40 


-69 19 57.5 


2.0 


X-rays (1) 


05 40 11.16 


0.15 


-69 19 57.79 


0.90 


Optical (2) 


05 40 10.980 


0.090 


-69 19 55.17 


0.50 


Optical (3) 


05 40 11.221 


0.090 


-69 19 54.98 


0.50 


X-rays (4) 


05 40 11.173 


0.120 


-69 19 54.41 


0.70 


Optical (5) 


05 40 11.160 


0.040 


-69 19 53.90 


0.20 


X-rays (6) 


05 40 11.202 


0.009 


-69 19 54.17 


0.05 


Optical (7) 



(1) Seward et al. (1984); (2) Manchester & Peterson (1989); (3) 
Caraveo et al. (1992); (4) Kaaret et al. (2001); (5) Serafimovich 
et al. (2004); (6) Livingstone et al. (2005); (7) this work 



4. Discussion 

4.1. The pulsar astrometry 

Unfortunately, the lack of a measurable proper motion for 
PSR B0540— 69 does not allow to search for possible con- 
nections between the pulsar kinematics, its polarisation 
properties, and the PWN morphology, which have been 
found, e.g. for the Crab and Vela pulsars (see Mignani 
2009c and references therein). The computed 3ct upper 
limit of 1 mas yr" 1 (see Sect. 3.1) on the PSR B0540-69 
proper motion sets a corresponding upper limit of ~ 250 
km s _1 on its transverse velocity. Although this value is 
lower than the peak of the pulsar transverse velocity distri- 
bution (e.g., Hobbs et al. 2005), it does not allow to claim 
a peculiarly low velocity for PSR B0540— 69 unless one 
lower the constraints on the measured upper limit to la. 
Indeed, several pulsars feature transverse velocities lower 
than 100 km s _1 including, e.g. the Vela pulsar for which 
a velocity of ~ 65 km s" 1 has been inferred from HST 
astrometry (Caraveo et al. 2001). 

The computed average coordinates of PSR B0540— 69, 
together with their associated errors, are listed in Tabled 
for comparison with those obtained from previous works. 
After its discovery by Einstein (Seward et al. 1984), the 
coordinates of PSR B0540— 69 were firstly revised by 
Manchester & Peterson (1989), from optical timing ob- 
servations, and later by Caraveo et al. (1992), from the 
astrometry of the optical counterpart identified in their 
NTT images. We note that the discovery of the pulsar 
in radio (Manchester et al. 1993) did not allow to ob- 
tain more precise coordinates. Indeed, its very low radio 
flux density at the LMC distance, S400 = 0.7 /iJy, makes 
it difficult to obtain an accurate radio timing position. 
Thus, the optical coordinates of Caraveo et al. (1992) were 
thereafter assumed as a reference. In particular, they were 
used to compute the pulsar X-ray and radio timing solu- 
tion by, e.g. Zhang et al. (2001), Cusumano et al. (2003) 



but, surprisingly enough, not by Finley et al. (1993) and 
Manchester et al. (1993). Chandra observations (Kaaret 
et al. 2001) yielded an updated X-ray position. The for- 
mer, and actually not the latter as the authors claim, was 
used as a reference by Johnston et al. (2004) for the pulsar 
X-ray/radio timing solution, and by de Plaa et al. (2003). 
A new optical position was obtained by Serafimovich et 
al. (2004) from HST data, while a new X-ray position was 
derived by Livingstone et al. (2005) based on Rossi-XTE 
timing observations. 

To visualise the difference and the improvement on the 
determination of the PSR B0540—69 position, we plotted 
all coordinates from Table[3]on the October 1995 WFPC2 
F555W image of the field (Fig. [I]). We see that the origi- 
nal pulsar coordinates from Seward et al. (1984) and from 
Manchester & Peterson (1989) fall « 3'/5 away from the 
pulsar position and, actually, out of the nebula shell. On 
the other hand, the coordinates of Caraveo et al. (1992) 
and Kaaret et al. (2001) are indeed within the nebula 
but they are clearly offset from the actual pulsar position, 
while those of Serafimovich et al. (2004) and Livingstone 
et al. (2005) are more consistent with the pulsar posi- 
tion, although the former have much larger errors. Thus, 
thanks to the sharp angular resolution of the WFPC2 and 
to the high astrometric accuracy of 2MASS, our coordi- 
nates are the most accurate obtained so far, providing a 
factor of ^ 4 improvement with respect to the most re- 
cently published values. In particular, our WFPC2 coor- 
dinates supersede both the optical coordinates of Caraveo 
et al. (1992), which were obtained using images with lower 
spatial resolution (O'.'13/pixel) and the GSC 1.0 (ct s - l'.'O; 
Lasker et al. 1990) as a reference catalogue, and those of 
Serafimovich et al. (2004), which were obtained using an 
early release of the GSC 2. 0, as well as the X-ray coordi- 
nates of Kaaret et al. (2001) and Livingstone et al. (2005), 
obtained through Chandra imaging and Rossi-XTE tim- 
ing observations, respectively. 
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Log(Frequency) (Hz) 



Fig. 1. 15" x 15" region of the PSR B0540-69 field 
observed with the HST/WFPC2 in the F555W filter 
(October 1995 observation). North to the top, East to 
the left. The crosses mark the different positions of the 
pulsar obtained from previous X-ray and optical observa- 
tions (black) and from our WFPC2 observation (white). 
The cross arms correspond to the la error in each coor- 
dinate. Different positions are numbered according to the 
reference publications (see Table ■ The pulsar is indi- 
cated by the arrow for clarity. 

The better determination of the PSR B0540-69 abso- 
lute position is important for follow-up non-imaging ob- 
servations. In the case of, e.g. spectroscopy, it will make 
it possible to accurately centre the pulsar in the slit when 
blind presets are made necessary, e.g. for ground-based 
observations due to the difficulties in resolving the pulsar 
through lower spatial resolution and seeing-dominated ac- 
quisition images. As a consequence, this will allow to use 
narrower slits without causing any loss of signal from the 
pulsar. This is important to minimise the background con- 
tamination from the supernova remnant which hampered 
ground-based spectroscopic observations performed so far 
(see, e.g. Serafimovich et al. 2004). 

4.2. The pulsar spectrum 

We used our multi-band WFPC2 photometry of PSR 
B0540— 69 to better investigate its optical SED. For con- 
sistency with previous works, we computed the interstel- 
lar extinction towards the pulsar using as a reference an 
E(B — V) = 0.2 and R — 3.1, which have been verified by 
several independent measurements (see Serafimovich et al. 
2004 and references therein) . We derived the extinction co- 
efficients in the WFPC2 filters from the extinction curves 
of Fitzpatrick (1999). The pulsar SED is shown in Fig.H 



Fig. 2. Optical spectral energy distribution of PSR 
B0540-69 derived from the available multi-band WFPC2 
photometry (Table Q}. Points are labelled according to the 
filter names. The dashed line is to the best fit power-law 
spectrum. 




10" 10 15 10 16 10" 10 ,a 10 19 



Frequency (Hz) 

Fig. 3. Optical spectral energy distribution of PSR 
B0540— 69 (points) compared with the power-law model 
(Kaaret et al. 2001) best-fitting the Chandra X-ray spec- 
trum (solid line) and its extrapolation in the optical do- 
main. Dashed lines correspond to a la uncertainty on the 
model parameters. 

after correction for the interstellar extinction. The plotted 
errors on the spectral fluxes are purely statistical and do 
not account for the systematic 0.02-0.04 magnitude error 
on the WFPC2 zero points (Heyer et al. 2004). A linear fit 
to the data points yields a spectral index qq = 0.70±0.07. 
As seen, our fit clearly confirms that there is no evidence 
for the "U-band excess" originally claimed by Middleditch 
et al. (1985). 
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We compared our best-fit power-law spectrum with 
that measured by Serafimovich et al. (2004). They de- 
rived a spectral index ao = 1.07io 19 which is somewhat 
steeper than that measured by us, although not formally 
incompatible. As we stated above, the power-law slope is 
only marginally affected by the difference in the assumed 
interstellar extinction correction. Indeed, applying the ex- 
tinction coefficients of Cardelli et al. (1989) to our spec- 
tral fluxes, as done by Serafimovich et al. (2004), yields 
a power law slope which is substantially identical to that 
obtained using the extinction coefficients of Fitzpatrick 
(1999). Thus, the difference in the power-law slope is most 
likely intrinsic to the photometry measurements. First of 
all, Serafimovich et al. (2004) fitted a lower number of 
points (see Table [lj and used the flux measured in the 
medium/narrow-band F547M and F658N filters, while we 
fitted a larger number of points and we only used the flux 
measured in wide-band filters. Moreover, while we used 
similar photometry apertures, we applied aperture cor- 
rections specific for the WFPC2, which are wavelength- 
dependent (Holtzman et al. 1995), and we applied the 
CTE correction which, apparently, was not accounted for 
by Serafimovich et al. (2004). 

The updated value of the PSR B0540— 69 optical spec- 
tral index is now more consistent with those of the other 
optically identified pulsars which all feature relatively 
steep spectral indexes (ao <; 0.4), while only the Crab 
and Vela pulsar feature a nearly flat power-law spectrum 
(see Mignani et al. 2007a for a comparison). In particu- 
lar, PSR B0540— 69 may be the pulsar with the steepest 
optical spectral index. We also note that the new value 
of ao = 0.70 ± 0.07 is very close to that of the X-ray 
spectral index ax — 0.83 ± 0.13 measured with Chandra 
(e.g., Kaaret et al. 2001). Interestingly, however, the op- 
tical spectrum is not the continuation of the X-ray one. 
Indeed, the HST spectral fluxes lie about a factor of 100 
below the low-energy extrapolation of the Chandra 0.1-10 
keV X-ray spectrum (see Fig. [3]). This difference can not 
be accounted for by possible problems in the subtraction 
of the PWN background or by other kinds of systematic 
effects, and confirms the presence of a spectral turnover 
between the optical and the X-ray bands (see also Fig. 14 
of Serafimovich et al. 2004). 

Spectral turnovers between these two energy bands are 
a common feature of the magnetospheric emission of many 
other rotation-powered pulsars, as indicated by the dif- 
ferences in their optical and X-ray spectral indexes (e.g., 
Mignani et al. 2007a). For comparison, we plotted in Fig. 
0] the optical, near ultraviolet (NUV), and near infrared 
(NIR) spectral fluxes for all rotation-powered pulsars with 
an optical counterpart, together with the model 0.1-10 keV 
X-ray spectra. In the case of the Crab and Vela pulsars, 
the spectral turnover is in the form of a single break in the 
optical-to- X-ray power-law spectrum. In the case of PSR 
B0540— 69, instead, a double break is required to join the 
optical and the X-ray SEDs, unless the actual interstellar 
extinction is larger by a factor of 2 with respect to current 
best estimates, which would make the spectrum flatten in 



the blue. This apparent double break in the optical-to- 
X-ray power-law spectrum might be also present in other, 
much closer, rotation-powered pulsars. Of course, for some 
pulsars (especially for the fainter ones) the comparison is 
hampered by the paucity of spectral flux measurements 
at NIR/optical/NUV wavelengths, which makes it diffi- 
cult to establish the presence of a power-law component. 
Bearing this caveat in mind, a possible double break can 
be recognised in the optical-to-X-ray power-law spectrum 
of the middle-aged pulsars PSR B0656+14 (wl00 kyrs) 
and, to a lesser extent, Geminga (« 340 kyrs), for both 
of which a power-law tail is hinted in the NIR. On the 
other hand, a single break is probably required for the old 
pulsars PSR B1929+10 (w 3 Myrs) and PSR B0950+08 
(« 18 Myrs), while only for the young PSR B1509-58 
the optical/NIR power-law seems to nicely fit the extrap- 
olation of the X-ray one. So, three different trends are 
recognised in our sample which spans four age decades. 
This means that the breaks in the power-law spectrum do 
not correlate with the neutron star age and, thus, are not 
indicative of any evolution in the emission processes in 
the neutron star magnetosphere. This is consistent with 
the lack of evidence for a power-law spectrum evolution, 
both in the optical (Mignani et al. 2007a) and in the X- 
rays (Becker 2009). The shape of the power-law optical- 
to-X-ray spectrum does not correlate with the dipole mag- 
netic field either with, e.g. the Crab and PSR B0950+08 
both featuring a single spectral break but having mag- 
netic fields different by two orders of magnitudes. It is 
possible that the number of observed spectral breaks is 
related to the geometry of the emission regions and to the 
particle distribution in the neutron star magnetosphere 
which, however, are difficult to reconstruct without hav- 
ing information on both the X-ray and optical light curve 
profiles. At the same time, is not possible to determine 
whether such double breaks are simply associated with a 
change in the power-law spectrum or they are associated 
with absorption processes in the neutron star magneto- 
sphere, a possibility proposed by Serafimovich et al. (2004) 
for PSR B0540-69. Observations in the far UV (FUV) 
would be fundamental to clarify this issue. Unfortunately, 
only few pulsars have been observed in this spectral re- 
gion by the EUVE satellite (Korpela & Bowyer 1998) but 
only Geminga (Bignami et al. 1996) and PSR B0656+14 
(Edelstcin et al. 2000) have been detected amongst the 
pulsars shown in Fig. 01 with their FUV fluxes being more 
or less compatible with the extrapolations of the black 
body components fitting the optical spectra. 

4.3. The pulsar polarisation 

Our value of the phase-averaged polarisation (PD = 
16% ± 4%) is consistent with the upper limit of 15% in- 
ferred by Middlcditch et al. (1987) from phase-resolved 
polarimetry of the pulsar. However, the measured PD is 
higher than PD ss 5% obtained from VLT phase-averaged 
polarisation observations by Wagner & Seifert (2000), for 
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which, however, no error estimate is given, so that its sig- 
nificance can not be assessed. Thus, ours is the only sig- 
nificant measurement (~ 4<t level) of the PSR B0540— 69 
phase-averaged polarisation obtained so far. 

Measurements of the phase-averaged polarisation de- 
gree have been obtained so far for all the young (Ss 10 000 
years old) rotation-powered pulsars with an optical coun- 
terpart, albeit with a different degree of confidence (see 
Mignani et al. 2007b and Slowikowska et al. 2009 for a 
summary). For the Crab, a value of PD = 9.8% ± 0.1% 
was inferred from phase-resolved polarimetry observations 
(Slowikowska et al. 2009), while image polarimetry obser- 
vations with the VLT yielded PD = 9.4% ± 4% for the 
Vela pulsar (Mignani et al. 2007b) and PD = 10.4% for 
PSR B1509-58 (Wagner & Seifert 2000), with the latter 
measurement being admittedly very uncertain and quoted 
with no error. The phase-averaged polarisation of PSR 
B0540-69, PD = 16% ± 4%, is thus consistent with the 
measurements obtained for the other young pulsars, as 
one might expect from their similar parameters, like spin 
down luminosity and/or magnetic field strength. Although 
the current data base is still extremely limited, with the 
measurement obtained for PSR B1509— 58 in wait for con- 
firmation and with the phase-resolved polarimetry obser- 
vations of the middle-aged (~ 100 kyrs) PSR B0656+14 
(Kern et al. 2003) covering only one third of the period, it 
nonetheless suggests that the phase-averaged optical po- 
larisation degree of rotation-powered pulsars is typically 
around 10%, with our value of PD for PSR B0540-69 
being possibly somewhat larger. 

On the theoretical side, the comparison of the observed 
phase-averaged polarisation degrees with the predictions 
of different pulsar magnetosphere models is complicated 
both by the degree of complexity of such models and by 
the limits in model simulations. However, for both the 
Crab and Vela pulsars, a detailed comparison with differ- 
ent pulsar magnetosphere models, showed that the mea- 
sured phase-averaged polarisation degree of the emerging 
optical radiation requires rather low degree of polarisation 
intrinsic to the source and/or strong depolarisation factors 
(Slowikowska et al. 2009; Mignani et al. 2007b). Recent 
3D outer- gap model calculations of optical and/or high- 
energy radiation with detailed incorporation of electron- 
positron gyration have been carried out so far for the Crab 
pulsar only (Takata et al. 2007). In general, the results 
seem to fit the measured optical polarization level (see 
Slowikowska et al. 2009) for some particular values of the 
viewing angle (the angle between the line of sight and the 
spin axis). Other models addressing the problem of po- 
larisation characteristics are relatively more simplified in 
terms of the model assumptions in use (Dyks et al. 2004; 
Petri & Kirk 2005). All models mentioned above - relevant 
for young energetic pulsars - rely on spatially extended 
sources of emission (either within the magnetosphere - 
Dyks et al 2004; Takata et al. 2007) or in the pulsar wind 
zone (Petri & Kirk 2005). As a consequence, the depolari- 
sation effects due to rotation, photon finite time of flight, 
and magnetic field structure are significant in all three 



models. This might explain the somewhat larger polarisa- 
tion degree of PSR B0540— 69. However, one should keep 
in mind that the shape of its optical light curve makes it 
more difficult to infer the geometry of the emission region 
with respect to the Crab and Vela pulsars, which provides 
crucial input parameters for model simulations. 

Interestingly, for both the Crab and Vela pulsars the 
polarisation position angle features a remarkable align- 
ment with the axis of symmetry of the X-ray structures 
(torus and jets) observed by Chandra, with the pulsar spin 
axis, and with the proper motion vector (Slowikowska et 
al. 2009; Mignani et al. 2007b). Unfortunately, for PSR 
B0540— 69 the scenario is less clear since no proper mo- 
tion has been measured so far and Chandra observations 
(Gotthelf & Wang 2000; Petre et al. 2007) could only par- 
tially resolve the morphology of the PWN, although its 
clear asymmetry might hint at the existence of a torus 
and, possibly, of a jet. A similar PWN morphology was 
also observed in the optical from HST/WFPC2 observa- 
tions (Caraveo et al. 2000). Most noticeably, the WFPC2 
images shows the existence of a bright emission knot south 
west of the pulsar which is apparently moving at a speed 
of ~ 0.04 c (De Luca et al. 2007) and which can be inter- 
preted as the head of a possible jet emitted by the pulsar. 
We note that the orientation of the pulsar polarisation vec- 
tor (PA = 22° ±12°) is consistent, accounting for possible 
projection effects, with that of the semi-major axis of the 
PWN (see Fig. [5]). Moreover, the pulsar phase-averaged 
polarisation vector is also possibly aligned with the ap- 
parent direction of motion of the knot (« 230° east from 
northjf]. We estimated the chance alignment probability 
between the two vectors to be ~ 0.04, small but still not 
negligible. If real, however, this peculiar alignment would 
indicate a possible physical connection between the pul- 
sar and the knot, as already proposed in De Luca et al. 
(2007), and would support the pulsar /jet scenario. The 
measurement of both the knot and of the PWN polarisa- 
tion structure (Mignani et al. in preparation) would fur- 
ther reinforce this scenario. Assuming for PSR B0540— 69 
the same PWN geometry as the Crab and Vela pulsars, 
the direction of the jet would suggest that a putative torus 
would indeed extend along the semi-minor axis of the neb- 
ula, and not along the semi-major one, as previously hy- 
pothesised from Chandra and HST observations (Gotthelf 
& Wang 2000; Caraveo et al. 2000). 



5. Summary 

By using the WFPC2 we performed a comprehensive 
study of the LMC pulsar PSR B0540-69 in the optical do- 
main, including astrometry, multi-band photometry, and 
polarimetry, down to the limits achievable with the pre- 
refurbished HST. 



We note that accounting for the proper motion of the LMC 
(Costa et al. 2009) and for the galactic rotation does not sig- 
nificantly affect this apparent alignment. 
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Fig. 5. Same as in Fig.Q]but zoomed over a 8" x 8" region. 
PSR B0540— 69 is labelled. The image has been smoothed 
over 3x3 pixel cells for a better visualisation of the ex- 
tended PWN emission. Isophotal contours (linear spacing) 
are drawn as solid lines. The thick lines show the measured 
la uncertainty (±12°) on the pulsar polarisation position 
angle, while the thin arrows show the assumed uncertainty 
(±30°) on the apparent proper motion direction of the 
bright emission knot in the PWN (De Luca et al. 2007). 

From the available multi-epoch images we obtained 
the tightest possible constraints on the PSR B0540— 69 
proper motion (< 1 mas yr _1 ) and transverse velocity 
(< 250 km s -1 ). Since the WFPC2 has been now de- 
commissioned and PSR B0540— 69 has not been observed 
yet with ACS, it is unlikely that our limit will be im- 
proved by HST before its retirement around 2015. Only 
the MICADO adaptive optics IR camera at the European 
Extremely Large Telescop^ (E-ELT), with an expected 
sensitivity limit of ~ 50/ias yr" 1 (Trippe et al. 2010), will 
be eventually able to measure the pulsar proper motion. 
From the same WFPC2 data set we measured the abso- 
lute position of PSR B0540— 69 with an accuracy of ~ 70 
mas, improving by a factor of 4 on the available X-ray 
coordinates. A further improvement could be achieved by 
HST using the refurbished ACS or the new Wide Field 
Camera 3 ( WFC3) , which have a larger field-of-view and 
a smaller pixel scale with respect to the WFPC2 mosaic, 
and the future Gaia astrometric catalogue which, how- 
ever, will not be released in its final version until 2019. 
From observations in six different filters we characterised 
the PSR B0540— 69 optical SED, as accurately as it could 
be done through broad-band imaging photometry with 
the HST, confirming that the pulsar spectrum is a mod- 
elled by a power-law (ao — 0.70 ±0.07) without evidence 



of a previously claimed spectral turnover in the U band. 
While multi-band photometry with ACS and WFC3 can 
certainly improve on our result, especially in the red and 
blue part of the spectrum, high-spatial resolution spec- 
troscopy with either of the two instruments or with the 
refurbished Space Telescope Imaging Spectrograph (STIS), 
is the only way to obtain a more detailed characterisation 
of the pulsar optical SED. Finally, from image polarime- 
try we derived a new measurement of the phase-averaged 
polarisation degree of PSR B0540-69 (PD = 16% ± 4%; 
PA = 22° ± 12°), with the polarisation vector apparently 
aligned with the PWN semi-major axis and with the ap- 
parent proper motion direction of its bright emission knot 
(De Luca et al. 2007). Multi-band polarimetry observa- 
tions with the refurbished, higher throughput, ACS will 
allow to obtain both a more precise measurement of the 
pulsar phase-averaged polarisation degree and of its po- 
sitional angle, and to investigate a possible dependence 
on the wavelength. Phase-resolved polarisation measure- 
ments in the optical, as well as in the X and 7-rays, are 
crucial for testing and constraining current models of pul- 
sar magnetosphere activity. Unfortunately, sensitive fast 
high-energy polarimeters, like eg. PoGOLite (Kamae et 
al. 2008), are still in the phase of developments, and they 
will not be available soon. For the time being the phase- 
resolved polarimetric studies in optical remain the only 
way to carry out such investigations. Therefore, phase- 
resolved optical polarimetry which, at present, can only 
be performed with guest instrument at ground-based tele- 
scopes, like Optima (Kanbach et al. 2009) and GASP 
(Collins et al. 2009), will be then unique and crucial to 
study the polarisation properties of noncoherent magne- 
tospheric radiation from PSR B0540— 69 and other op- 
tical pulsars. In particular, given the faintness of PSR 
B0540— 69 in the radio band, polarisation measurements 
in the optical domain are the only way to carry out such 
investigations. 

About 25 years after its discovery, PSR B0540— 69 is 
still the only extra galactic optical pulsar identified so far. 
Thus, it is a unique target for optical observations with 
the present ground and space-based facilities. More extra 
galactic neutron stars, however, are expected to be dis- 
covered in the optical, both in the Magellanic Clouds and 
beyond, after the advent of the E-ELT which holds the 
potential of detecting Crab-like optical pulsars up to the 
distance of M31 (Shearer 2008). 
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Fig. 4. Same as Fig. [3] but for all rotation-powered pulsars with an optical counterpart and flux measurements in at 
least two bands. Optical flux values are taken from the compilation in Mignani et al. (2007a). X-ray spectral models 
are taken from the source publications: Crab (Willingale et al. 2001), PSR B1509— 58 (Gaensler et al. 2002), Vela 
(Manzali et al. 2007), PSR B0656+14 (De Luca et al. 2005), Geminga (Caraveo et al. 2004), PSR B1929+10 (Becker 
et al. 2006), PSR B0950+08 (Becker et al. 2004). Left panels: young pulsars (r < 10 000 years). Right panels: middle- 
aged and old pulsars (r > 100 000 years). In each column objects are sorted according to their spin-down age and 
labelled. Different X-ray spectral components are shown by the dotted red, green (black body), and blue (power-law) 
lines while the solid blue lines show the composite spectra. Only best-fits are plotted for clarity. 



